Ultrafast carrier and phonon dynamics in BioSe^ crystals 
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Ultrafast time-resolved differential reflectivity of Bi 2 Se3 crystals is studied using op- 
tical pump-probe spectroscopy. Three distinct relaxation processes are found to con- 
tribute to the initial transient reflectivity changes. The deduced relaxation timescale 
and the sign of the reflectivity change suggest that electron-phonon interactions and 
defect-induced charge trapping are the underlying mechanisms for the three pro- 
cesses. After the crystal is exposed to air, the relative strength of these processes is 
altered and becomes strongly dependent on the excitation photon energy. 
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Recently there has been an emerging interest in Bi-based compounds, such as Bii_ z Sb z , 
Bi2Se3, and Bi2Te3, as exemplars of a new class of quantum matter known as a topologi- 
cal insulator-. Among these compound materials, Bi 2 Se3 stands out owing to its relatively 
simple topological surface states and its sizable bulk band gap of ~0.3 eV- 1 ^, which hold 
promise for applications in future spintronic and quantum computational devices. Indeed, 
in the past year the electronic properties of Bi2Se3 have attracted considerable attention^~- 
as field effect transistor type devices have been demonstrated experimentally- - — . However, 
in order to implement room-temperature operation of Bi2Se3-based devices, it is essential 
to understand the dynamics of phonons and excited carriers in this material^, particularly 
the electron-electron, electron-phonon and phonon-phonon interactions. In this Letter, we 
report the ultrafast time-resolved optical spectroscopy study of Bi 2 Se3 crystals in both the 
time domain and the energy domain. Our measurements reveal three underlying relaxation 
processes in the transient response of Bi2Se3, each associated with different physical mech- 
anisms. It is also shown that the relative strength of these processes is sensitive to air 
exposure of the samples. The observed charge trapping and air doping effects are likely due 
to the presence of Se vacancies, a major issue material scientists working to use the novel 
properties of Bi 2 Se3 will face in the near term. 

The Bi 2 Se3 single crystals studied in this work were synthesized via the Bridgman method 
at Purdue University and Fudan University. During crystal growth, the mixture of high 
purity elements was first deoxidized and purified by multiple vacuum distillations, and then 
heated to 850-900 °C for 15 hours, followed by a slow cool down under a controlled pressure 
of Se to compensate for possible Se vacancies. Afterwards, the samples were zone refined 
at a speed of 0.5-1.5 mm/hour with a linear temperature gradient set to 4-5 °C/cm, until 
a temperature of 670 °C was reached. The as-grown Bi 2 Se3 crystals from both groups are 
naturally n-doped due to remnant Se vacancies^. Hall measurements show typical carrier 
densities over 1 x 10 18 cm -3 at ambient conditions. 

In the pump-probe experiment, the transient reflectivity changes, AR/R, are measured 
at room temperature by employing a Ti:Sapphire laser with a repetition rate of 76 MHz. 
The laser is capable of producing ~150 fs-wide pulses in a wavelength range from ~720 nm 
(1.72 eV) to ~940 nm (1.32 eV). Cross-polarized (linear) pump and probe beams of the 
same wavelength are focused onto the sample in a spot ~100 /im in diameter. To avoid 
heating damage, the typical pump light fluence used in this work is 2.5 /iJ/cm 2 . 
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FIG. 1. Time-resolved response of Bi2Se3 at room temperature for the pump-probe photon en- 
ergy 1.63 eV over long (main panel) and short (inset) timescales. A schematic drawing of the 
measurement configuration is also shown in the inset. 

Figure [T] shows the measured AR/R with the pump-probe wavelength centered at 760 
nm (1.63 eV) over long (main panel) and short (inset) timescales. As one can see, AR/R 
undergoes a sharp rise within ~1 ps, followed by a decay. In addition, two apparent oscilla- 
tory behaviors are superimposed on the AR/R profile: One appears within the initial ~15 
ps with a higher frequency as evidenced in the small periodic "wiggles" in the inset to Fig. [TJ 
while the other arrives after ~6 ps and decays in a few periods with a much lower frequency 
but a larger amplitude (see the main panel of Fig. [T]). The low frequency oscillation has 
been widely observed in ultrafast time-resolved optical measurements of other systems^ and 
has been attributed to coherent longitudinal acoustic phonons generated by the pump light. 
To study the high frequency oscillation, we first extract it by subtracting the non-oscillatory 
component from the differential reflectivity data, and then fit it as an exponentially damped 
oscillation^. As shown in the inset to Fig. El our data can be fit well by this method 
with a single frequency v centered at 2.13 THz and a damping rate of T ~ 0.33 THz. Here 
the deduced v value is consistent with the energy of the A\ g longitudinal optical phonon in 
Bi 2 Se3 obtained by Raman spectroscopy^, while the damping rate represents the timescale 
of the phonon-phonon interactions. Moreover, v and T are found to be independent of the 
excitation photon energy E p h within our experimental error. 

Now let us discuss the non-oscillatory transient response of Bi 2 Se 3 in the first few picosec- 
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FIG. 2. (color online) AR/R signal of Bi2Se3 as a function of delay time at E p h = 1.63 eV. The best 
fit to data (red solid line) is obtained by using a model consisting of three independent exponential 
relaxation processes with time constant t± (short dashed line), r 2 (dashed-dot line), and T3 (dashed 
line). Inset: Extracted oscillations from the measured AR/R and the best fit to data (red solid 
line) as an exponentially damped oscillation. 

onds after pumping. Figure [2] plots AR/R as a function of delay time measured at room 
temperature and E p h = 1.63 eV. We note that the reflectance rise time in Bi 2 Se3 (~1 ps) 
is much longer than that expected for typical semiconductors and metals (~300 fs)^. This 
timescale corresponds to the delayed thermalization of the carriers via electron-electron in- 
teractions. Such a slow rise may indicate the presence of a relaxation process with a negative 
amplitude (^2), as illustrated by the dashed-dot line in Fig. [2j In addition, as further evi- 
denced in Fig. [31 two positive processes also contribute to the transient reflectivity changes, 
with distinctive time constant n and r 3 . Thus, to obtain the best fit to the data in Figs. [2] 
and [3] (the red solid lines), we employ a model consisting of three independent exponential 
relaxation processes^. The deduced time constants, T\ = 0.11 ± 0.02 ps, r 2 = 0.32 ± 0.03 
ps, and T3 = 2.3 ± 0.2 ps, are independent of the pump intensity and the excitation photon 
energy available with our experimental set-up. 

What is the possible origin of the three relaxation processes? The first process is charac- 
terized by a fast relaxation (n = 110 fs), as the excited carriers lose their high excess energy 
rapidly through the electron-phonon (particularly electron-optical-phonon) interactions^ - — . 
The third process, on the other hand, represents a relatively slow relaxation (r 3 = 2.3 ps), 



4 



in 

o 
5, 
CH 

< 

6 
4 
2 

-2 

1 2 3 1 2 3 

Delay time (ps) 

FIG. 3. (color online) Side-by-side comparison of the AR/R signals of Bi2Se3 exposed to air and 
in vacuum at three different excitation photon energies. Traces are offset for clarity. The red solid 
lines are the best fits to data with three independent exponential relaxation processes. 

which might be attributed to electron-acoustic-phonon interactions^ 1 ^. The co-existence of 
a sub-ps and a few-ps electron-phonon relaxation processes due to different phonon modes 
has been reported previously, and quantitative agreement is achieved in the extracted relax- 
ation timescales^ 1 ^. The second process with a negative amplitude A2, however, is more 
subtle. Such a negative-amplitude process was previously observed in strongly correlated 
electron systems^, as well as in disordered low-temperature molecular-beam-epitaxy grown 
III-V semiconductors^ 1 ^. One possibility is that this process is associated with the ultra- 
fast trapping of electrons by Se vacancies. This trapping mechanism could cause a negative 
change in reflectivity and limit the relaxation time to less than 1 psM-, in agreement with 
our observations. 

An important material issue for Bi2Se3 crystals is unintentional charge doping from air, 
which could modify or even dominate its electronic properties^. To pursue this issue, we 
perform comparative studies of the time-resolved response of Bi 2 Se3 exposed to air and in 
vacuum at different excitation photon energies. To prepare an "in-vacuum" sample, we 
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FIG. 4. | ^4.2 / ^4.1 1 and A^j A\ as a function of E p h for Bi2Se3 crystals exposed to air (filled dots) and 
in vacuum (open squares). Solid and dashed lines are the guides to the eye. 



first rapidly cleave it to obtain a mirror-like surface and load it into a sealed measurement 
chamber within 20 seconds, the chamber is then pump down to ~10~ 5 mbar in ~20 mins. 
As shown in Fig. [31 a clear difference in AR/R is evident as a function of E p h- Since 71,2,3 
exhibit almost no change when the samples are exposed to air, the observed differences are 
mainly from the relative changes in the amplitudes ^1,2,3- Figure H] summarizes the E p h 
dependence of | and A3/A1 for the measurements in air and in vacuum. While the 

relative strength of the three relaxation processes in vacuum is only weakly dependent on 
E p h, an abrupt change occurs for the samples exposed to air when E p h ~ 1.47 eV. Air 
doping may raise the Fermi level in Bi 2 Se3 and affect the band bending near the crystal 
surface. This alters the spatial electron-hole distribution function^ 2 ., which impacts the 
ultrafast dynamics of the material. Further experimental and theoretical work is needed to 
understand our observation. Finally, we find that the extracted A\ g optical phonon energy 
(v = 2.13 THz) from the in-vacuum samples is slightly, but consistently, lower than that 
obtained for the samples exposed to air (2.15 THz) at all E p h measured. Given that the 
theoretical value of this energy is 2.03 THz for a perfect Bi 2 Se 3 crystal^, Se vacancies may 
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also underlie the observed difference. 

In conclusion, we have performed ultrafast time-resolved optical spectroscopy measure- 
ments on Bi2Se3 crystals and observed the contributions from three exponential relaxation 
processes with distinctive time constants. The origin of the three processes can be attributed 
to electron-phonon interactions and defect-induced charge trapping. Our work also shows 
that air exposure can significantly effect the carrier and phonon dynamics in Bi2Se3 crystals 
by altering the relative strength of these three processes. 
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REFERENCES 

1 Fot a recent review, see X.-L. Qi and S.-C. Zhang. larXiv: 10 08. 2026, 
2 Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H. Lin, A. Bansil, D. Grauer, Y. S. Hor, R. 
J. Cava, and M. Z. Hasan, Nature Physics 5, 398 (2009). 

3 H. Zhang, C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and S.-C. Zhang, Nature Physics 5, 438 
(2009). 

4 Y. S. Hor, A. Richardella, P. Roushan, Y. Xia, J. G. Checkelsky, A. Yazdani, M. Z. Hasan, 
N. P. Ong, and R. J. Cava, Phys. Rev. B 79, 195208 (2009). 

5 A. D. LaForge, A. Frenzel, B. C. Pursley, T. Lin, X. Liu, J. Shi, and D. N. Basov, Phys. 
Rev. B 81, 125120 (2010). 

6 N. P. Butch, K. Kirshenbaum, P. Syers, A. B. Sushkov, G. S. Jenkins, H. D. Drew, and J. 
Paglione, Phys. Rev. B 81, 241301(R) (2010). 

7 J. G. Analytis, R. D. McDonald, S. C. Riggs, J.-H. Chu, G. S. Boebinger, and I. R. Fisher, 
larXiv:1003T7T3l 

8 J. Chen, H. J. Qin, F. Yang, J. Liu, T. Guan, F. M. Qu, G. H. Zhang, J. R. Shi, X. C. 
Xie, C. L. Yang, K. H. Wu, Y. Q. Li, and L. Lu, larXiv: 1003. 1531 

9 H. Steinberg, D. R. Gardner, Y. S. Lee, and P. Jarillo-Herrero, larXiv: 1003.31371 
10 J. G. Checkelsky, Y. S. Hor, R. J. Cava, and N. P. Ong, larXiv: 1003.38831 
11 J. Shah, Ultrafast spectroscopy of semiconductors and semiconductor nano structures 

7 



(Springer, 1999). 

J. Qi, J. A. Yan, H. Park, A. Steigerwald, Y. Xu, S. N. Gilbert, X. Liu, J. K. Furdyna, S. 
T. Pantelides, and N. Tolk, Phys. Rev. B 81, 115208 (2010); and the references therein. 
A. Q. Wu, X. Xu, and R. Venkatasubramanian, Appl. Phys. Lett. 92, 011108 (2008). 
W. Richter, H. Kohler, and C. R. Becker, Physica Status Solidi B-Basic Research 84, 619 
(1977). 

S. Rast, M. L. Schneider, M. Onellion, X. H. Zeng, W. Si, X. X. Xi, M. Abrecht, D. Ariosa, 
D. Pavuna, Y. H. Ren, G. Liipke, and I. Perakis, Phys. Rev. B 64, 214505 (2001). 
Y. Xu, M. Khafizov, L. Satrapinsky, P. Kus, A. Plecenik, and R. Sobolewski, Phys. Rev. 
Lett. 91, 197004 (2003). 

P. B. Allen, Phys. Rev. Lett. 59, 1460 (1987). 

T. Taneda, G. P. Pepe, L. Parlato, A. A. Golubov, and R. Sobolewski, Phys. Rev. B 75, 
174507 (2007). 

J. Demsar, V. K. Thorsm0lle, J. L. Sarrao, and A. J. Taylor, Phys. Rev. Lett. 96, 037401 
(2006). 

S. Gupta, M. Y. Frankel, J. A. Valdmanis, J. F. Whitaker, G. A. Mourou, F. W. Smith, 
and A. R. Calawa, Appl. Phys. Lett. 59, 3276 (1991). 

J. G. Wang, C. J. Sun, Y. Hashimoto, J. Kono, G. A. Khodaparast, L. Cywihski, L. J. 
Sham, G. D. Sanders, C. J. Stanton, and H. Munekata, J. Phys: Condens. Matter 18 R501 
(2006). 

T. Dekorsy, G. C. Cho, and H. Kurz, Chapter 4 in Light Scattering in Solids VIII: 
Fullerenes, Semiconductor Surfaces, Coherent Phonons (Topics in Applied Physics, Vol- 
ume 76) (Springer, 2000). 



S 



